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We investigate the electronic structures of iron-based superconductors having perovskite-like 
blocking layers, Ca4Al20eFe2P2 and Ca4Al20eFe2As2, from first principles. Ca4Al20eFe2P2 
is found to have two hole- like Fermi surfaces around T, and one hole- like Fermi surface around 
M in the unfolded Brillouin zone. This is in contrast with LaFePO, where no Fermi surface is 
found around M. The relationship between their band structures and measured transition tem- 
peratures of superconductivity is discussed. The number of Fermi surfaces in Ca4Al20eFe2P2 is 
also different from that in Ca4Al20eFe2As2, in which only one Fermi surface is formed around 
r. Analysis using maximally localized Wannier functions clarifies that the differences between 
their band structures originate mainly from the difference in pnictogen height. We then ana- 
lyze the alloying effect on the electronic structure of Ca4Al2 0eFe2AsP. It is found that this 
electronic structure is similar to those of Ca4Al20eFe2P2 and Ca4Al20eFe2As2 with average 
crystal structures, though Ca4Al20eFe2AsP contains pnictogen height disorder. We calculate 
the generalized susceptibility for Ca4Al206Fe2(Asi_ cc P cc )2 and clarify the factors determining 
its tendency. 

KEYWORDS: Fe-based superconductor, Ca4Al20eFe2P2, Ca4Al206Fe2As2, LaFePO, electronic structure, 
first-principles calculation, maximally localized Wannier function 



1. Introduction 

Since the discovery of iron-based superconductors in 
LaFePO with its transition temperature T c of 5 K 1 ) and 
in LaFeAsO with T c = 26 K, 2 ) many attempts to clar- 
ify the material properties and to realize higher T c have 
been continued, and Fe-based superconductors with var- 
ious chemical formulae and crystal structures have been 
reported, among which SmFeAsO has shown the highest 
T c of 55 K to date is realized. 3 ) 

An Fe-based superconductor comprises two- 
dimensional iron (Fe)-anion (A) layers. The elec- 
tronic bands near the Fermi level have a strong Fe 
d character, 4 7 ) which should be responsible for the 
superconductivity. Lee et al 8 ^ plotted the T c of various 
Fe-based superconductors as a function of the A-Fe-A 
bond angle a and demonstrated that the T c shows a 
peak at nearly 109.47°, for which Fe^4 forms a regular 
tetrahedron. Much attention has thus been paid to the 
strong correlation between T c and the local geometry 
of an FeA^ tetrahedron, namely, the lattice constant 
a, the Fe-A bond length d, the bond angle a, and the 
anion height Ha- Two of these four geometry constants 
are independent. A small a or d leads to a small density 
of states at the Fermi level and is thus expected to be 
unfavorable for superconductivity. On the other hand, 
a large Ha allows a robust hole Fermi surface around 
(7r, 7r), 9 ) which interacts with electron Fermi surfaces and 
is thus expected to be favorable for superconductivity. 10 ) 
The measured T c 's, however, are not necessarily higher 
at larger ft^'s. Mizuguchi and coworkers 11 ' 12 ) provided a 
plot of T c as a function of ft ,4, which obeys a symmetric 
curve with a peak at approximately 1.38 A, common to 



1111-, 122-, 111-, and 11-type superconductors. 

Since the reports for Sr4Sc2 06Fe2P2 13 ^ and 
Sr4V2 0eFe2As2, 14 ) Fe-based materials with perovskite- 
type blocking layers have been known to show super- 
conductivity for various combinations of their formulae 
and thickness. 15 23) They include Ca 4 Al 2 6 Fe2As2 
(a = 3.713 A and c = 15.407 A, a = 102.13°, and 
ft As = 1.500 A) with T c = 28.3 K 20 ) and Ca4Al 2 6 Fe2P2 
(a = 3.692 A and c = 14.934 A, a = 109.45°, and 
ftp = 1.306 A) with T c = 17.1 K. 2 °) 

Among the Fe-based superconductors reported so far, 
Ca4Al206Fe2As2 and Ca4Al20eFe2P2 are of particular 
interest since they have rather small a's, owing to the 
small ion radius of Al 3+ , and Ca4Al2 0eFe2As2 has a 
rather large ftAs- Furthermore, it is important to an- 
swer why the T c of Ca4Al2 0eFe2P2 is higher than that 
of LaFePO to elucidate the mechanism of superconduc- 
tivity in systems consisting of FeP layers. Our previ- 
ous work 24 ) revealed that a has a strong impact on 
the band rearrangement and the Fermi surface topology 
in Ca4Al2 0eFe2As2. Subsequently, Usui and Kuroki 25 ) 
pointed out that the superconductivity of an Fe-based 
superconductor is optimized for a of a regular tetrahe- 
dron, which maximizes the hole Fermi surface multiplic- 
ity. 

The superconducting properties of Fe-based systems 
depend sensitively on multiple factors related to the crys- 
tal structure and formula, as outlined above. It is thus 
important to investigate the electronic properties of such 
systems with the crystal structure and/or formula varied, 
focusing on the behavior of Fermi surfaces and electronic 
orbitals in the vicinity of the Fermi level. In the present 
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study, we therefore perform first-principles electronic 
structure calculations for Ca4Al206Fe2(Asi_ :E P :E )2. We 
first compare the band structures of Ca4Al2 0eFe2P2 
and Ca4Al2 06Fe2As2, and determine the origins of the 
differences in their band structures by comparing the 
transfer integrals between localized electronic orbitals. 
We extract the transfer integrals by constructing max- 
imally localized Wannier functions, 26 ) which provide a 
picture of localized electronic orbitals in solid. The re- 
sults of this analysis suggest that the band structures 
of Ca4Al206Fe2(Asi_ iC P x )2 are determined mainly by 
crystal structure, not by chemical composition. We then 
demonstrate that electronic structure calculation for 
Ca4Al20eFe2AsP corroborates this idea by unfolding its 
band structure. We calculate the generalized susceptibil- 
ity xo? which correlates with superconductivity in gen- 
eral, of Ca4Al206Fe2(Asi_ :E P :c )2 and clarify the factors 
determining the behavior of xo- 

2. Computational Details 

We used the computational code QMAS 27 ) based on 
the projector augmented- wave method, 28 ) which has 
been applied to the study of the ground state properties 
of LaFeAsO, 6,29 ) SrFe 2 As 2 , 30) and Ca 4 Al20 6 Fe2As2. 24) 
The PBE exchange-correlation functional 31 ) within the 
generalized gradient approximation (GGA) was adopted. 
The pseudo wave functions were expanded in plane waves 
with an energy cutoff of 40 Ry. We employed a 12x12x4 
k point mesh for systems with conventional tetragonal 
cells and an 8 x 8 x 4 k point mesh with super cells. The 
electronic band structure in the vicinity of the Fermi level 
is analyzed in detail using the maximally localized Wan- 
nier function technique. 26 ) 

For the Ca 4 Al2O 6 Fe2P2 20 ' 32) and Ca 4 Al2O6Fe 2 As2 20) 
crystal structures, the experimental values of tetragonal 
lattice constants and atomic coordinates are used in the 
present work. 

In the analyses of the alloying effect on 
Ca4Al20eFe2AsP, we unfold its band dispersion 
obtained for a supercell, following the procedures 
proposed by Ku et al 33 ^ The unfolded band dispersion 
is given as a sum of weighted delta functions: 

A kn , kn (u) = \(kn\KJ)\ 2 S(uj - e K j): (1) 
j 

where K is the wave vector in the folded Brillouin zone 
corresponding to the wave vector k in the unfolded Bril- 
louin zone. J is the band index for the supercell and n 
denotes the Wannier function in the primitive cell. \KJ) 
is the eigenstate of the Kohn-Sham Hamiltonian for the 
supercell, whereas \kn) is the Fourier transform of the 
Wannier function in the primitive cell. We are able to 
unfold 20 Fe d bands to 5 bands by choosing the proper 
gauge of the Wannier functions. 7 ) 

We calculate the generalized susceptibility xo of 
Ca4Al20eFe2P2 and Ca4Al20eFe2As2 with gradually 
varying their crystal structures. Its expression for a sys- 



tem with time reversal symmetry is given by 

/ \ 1 \ ^ frnk ~ fnk+q / \ 

xo(q) = o 2^ ; rr^' ( 2 ) 

m,n,k 

where / m & is the occupation number and e m k is the en- 
ergy eigenvalue, xo reflects the nesting property of the 
system and is expected to correlate with superconduc- 
tivity. We evaluate eq. (2) as a summation for 40 x 40 x 4 
k points by interpolating the band dispersion using the 
Wannier functions of Fe d bands. 

3. Results and Discussion 

3.1 Electronic properties of Ca/±Al20§Fe2P2 in compar- 
ison with those of Ca±Al20§Fe2As2 

The nonmagnetic electronic band structures of 
Ca4Al20eFe2P2 and Ca4Al20eFe2As2 obtained with the 
first-principles calculation are shown as circles in Fig. 
1(a). In each of the systems, the unit cell contains two 
Fe atoms and ten bands having a strong Fe d charac- 
ter lie near the Fermi level. The most prominent differ- 
ence in the calculated band structure between these two 
systems is the multiplicity of Fermi surfaces around T: 
Ca4Al20eFe2P2 has three, whereas Ca4Al2 0eFe2As2 has 
two. We note here that LaFePO has two Fermi surfaces 
around T, which would explain why Ca4Al20eFe2P2 
has a higher T c than LaFePO. Figure 1 (b) shows the 
Fermi surfaces of Ca4Al 2 06Fe2P2 in an undoped case 
and doped cases with ±0.1 electrons per Fe atom within 
the rigid band approximation. Another noticeable differ- 
ence is seen in the height of pnictogen (Pn) p bands at 
M. The top of P p bands in Ca4Al2 0eFe2P2 lies close to 
— 1.3 eV at M, whereas that of As p bands is not at M 
and those bands are below —1.9 eV. 

To analyze the differences between the two com- 
pounds, the Wannier function technique is useful, since 
it enables us to evaluate the effective parameters of low- 
energy electronic states. 34 ) We constructed maximally lo- 
calized Wannier functions from the ten Fe d bands (d 
model) and plotted the interpolated band dispersion as 
curves in Fig. 1(a). It is found that the Wannier func- 
tions accurately reproduce the DFT-GGA band struc- 
tures near the Fermi level. The Wannier function of the 
Fe d X 2-Y 2 (d X y) orbital in Ca4Al206Fe2P2 is shown in 
Fig. 1(c). [The X- and F-axes refer to the conventional 
unit cell containing two Fe atoms, whereas the rotation 
of them by 45° defines the x- and ?/-axes. (see Fig. 2)] 
Although the shape of this Wannier function is mainly 
of the dx^-Y 2 character, it involves contributions from 
P p orbitals. We observed that the positive (blue) lobe of 
the dx^-Y 2 orbital is pushed down the c direction to hy- 
bridize with the P p orbitals below the Fe layer, while the 
negative (yellow) lobe is pushed up to hybridize with the 
P p orbitals above, as shown in Fig. 1(c). We confirmed 
that the shape of the corresponding Wannier function in 
Ca4Al20eFe2As2 (not shown) has smaller contributions 
from the As p orbitals. This reflects (i) the Fe-P dis- 
tance, which is smaller than the Fe-As distance, and (ii) 
the energy separation of the Fe d and the As p bands in 
Ca4Al20eFe2As2, which is larger than that of the Fe d 
and the P p bands in Ca4Al2 0eFe2P2- 
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Fig. 2. (Color online) Schematic illustration of transfer integrals 
of a dpp model. Solid lines represent unit cells. Fe d X 2_ Y z or - 
bitals and Pn px and pz orbitals are shown. We define the 
x- and y-axes by rotating the original X- and F-axes by 45°. 
Transfer integrals are shown as red arrows. ti(i = 1, 2, 3) is that 
between the d X 2_ Y z orbitals at the z-th nearest neighbors, while 
tj (j = X, Z) is that between the pj orbital and its nearest- 



neighboring d X 2_ Y 2 orbital. 



Fig. 1. (Color) (a) Nonmagnetic electronic band structures of 
Ca4Al 2 6 Fe2P2 and Ca4Al 2 06Fe 2 As2. Circles represent DFT- 
GGA bands and curves represent bands interpolated using the 
maximally localized Wannier functions of d models. The origins 
of energy are set to the respective Fermi levels, (b) Fermi sur- 
faces for hole-doped (left), undoped (middle), and electron-doped 
(right) Ca4Al2C>6Fe2P2- (c) A bird's-eye view of maximally lo- 
calized Wannier function of Fe d X 2_ Y 2 orbital is shown on the 
left. The Wannier function viewed along the c axis is shown on 
the right. The a and b directions are parallel to the X- and Y- 
axes, respectively. Brown and purple balls represent Fe and P 
atoms, respectively. These figures were drawn with VESTA. 35 ) 



To simplify the analysis of the electronic band struc- 
ture, we unfolded the interpolated 10-band dispersion of 
d models following the procedure proposed by Kuroki 
et al 7 ^ The resultant 5-band structures obtained with 
the lattice parameters and the atomic coordinates fixed, 
while varying the chemical composition of the system, 
are shown in Fig. 3. It is found that the overall band 
structures with the same crystal structure look quite 
similar, especially in the vicinity of the Fermi level. For 
the Ca4Ai206Fe2P2 (Ca4Al2 0eFe2As2) crystal struc- 
ture, there exist three (two) Fermi surfaces around T 
in the original Brillouin zone irrespective of chemical 
composition, which correspond to the two (one) a hole 
pockets around (0,0) and one 7 hole pocket around 
(7r,7r), as shown in Fig. 3. The numbers of Fermi sur- 
faces around (0, 0) and (7r, it) in the unfolded Brillouin 
zone for Ca4Al20eFe2 As2 and Ca4Al20eFe2P2 with their 
own crystal structures are summarized in Table I. We 
observe in Fig. 3 that the calculated band widths of d 



models for Ca4Ai2 06Fe2As2 with the Ca4Ai2 06Fe2As2 
and Ca4Al2 06Fe2P2 crystal structures are 4.64 and 5.68 
eV, respectively. Those for Ca4Al2 06Fe2P2 with the 
Ca4Al206Fe2As2 and Ca4Al20eFe2P2 crystal structures 
are 4.29 and 5.28 eV, respectively. The Ca4Al 2 6 Fe2P2 
crystal structure gave a larger band width than the 
Ca4Al206Fe2As2 crystal structure when the same chem- 
ical composition was used, since the former has a 
smaller a and a smaller hp n . On the other hand, the 
Ca4Al20eFe2P2 composition gave a smaller band width 
than the Ca4Al20eFe2As2 composition when the same 
crystal structure was used, since the ion radius of the P 
atom is smaller than that of the As atom. 



Table I. Numbers of Fermi surfaces around (0, 0) and (7r, 7r) in un- 
folded Brillouin zone for Ca4Al2C>6Fe2As2 and Ca4Al20eFe2P2- 
Those for LaFeAsO and LaFePO are also shown for comparison. 
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It was demonstrated in our previous work 24 ) that the 
Fe d x z-Y 2 (dxy) orbital plays an important role in the 
band rearrangement of Ca4Al20eFe2As2 when the lo- 
cal geometry of the FePn4 tetrahedron is varied. For 
the quantitative analysis of the electronic structure of 
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Fig. 3. (Color online) Unfolded 5-band structures of d mod- 
els for Ca4Al2C>6Fe2P2 and Ca4Al20eFe2As2 as solid and 
dashed curves, respectively. The upper (lower) panel shows 
the band structures obtained using the crystal structure of 
Ca4Al2C>6Fe2P2 (Ca4Al20eFe2As2). The origins of energy are 
set to the respective Fermi levels. 

Ca4Al20eFe2P2 along this line, we constructed dpp mod- 
els of the four systems mentioned above and examined 
how their chemical composition and crystal structures 
affect the electronic band structures. Each of the dpp 
models consists of 10 Fe d, 6 Pn p, and 18 O p or- 
bitals. The transfer integrals and orbital energies of the 
dx 2 -Y 2 orbitals, which are equivalent to the xy orbitals, 
and the Pn p orbitals are shown in Table II. Although 
|ti| ^> |^2 1 in the dpp models, t\ in the d models are 
small in magnitude; consequently, t\ and t<i are compa- 
rable in the d models. Although the Wannier functions 
in the d model are spatially extended and have con- 
siderable weight at the pnictogen site, the d Wannier 
functions in the dpp model are localized and atomic- 
orbital-like. 36 ' 37 ) Hence, we can say that t\ in the dpp 
model contains only direct hopping between the neigh- 
boring Fe sites, whereas t\ in the d model contains in- 
direct hopping via the pnictogen site as well. 24 ) This in- 
dicates that the indirect electron hopping of Fe-Pn-Fe 
partly cancels the direct hopping between the nearest- 
neighboring d X 2-Y 2 orbitals. In the comparison of dpp t's 
for the same chemical composition but different crystal 
structures, the difference between ti(i = l,2)'s is small 
compared with their amplitudes. On the other hand, tx 
and tz are much larger in the Ca4Al2 0eFe2P2 crystal 
structure than in the Ca4Ai2 06Fe2As2 crystal structure 
for both chemical compositions. This is because the hy- 
bridization between the dx 2 -Y 2 an d p orbitals is stronger 
for a smaller hp n . Since tx and tz in the dpp models can 



be regarded to be incorporated into t's in the d models, 
the differences in band dispersion between the systems 
with the same chemical composition but different crystal 
structures come from the differences between the indi- 
rect hopping of Fe-Pn-Fe. The two-dimensional shape of 
the dx 2 -Y 2 orbital allows a larger hybridization with the 
pz orbital than with the px orbital, leading to tz much 
larger than tx for all the four systems. In the compar- 
ison of dpp Vs for the same crystal structure but dif- 
ferent chemical compositions, those for Ca4Al2 0eFe2As2 
are basically larger than those for Ca4Al20eFe2P2 owing 
to the ion radius of the As atom larger than that of the P 
atom, tx is, however, larger in Ca4Al2 06Fe2P2 than in 
Ca4Al20eFe2As2 if their crystal structures are the same. 
This is probably because the p electrons of a Pn atom 
are depleted for the hybridization between the pz orbital 
and its neighboring dx 2 -Y 2 orbitals. 

As clarified above, the differences between the band 
structures of Ca4Al 2 06Fe2P2 and Ca4Al 2 06Fe2As2 with 
their own crystal structures, shown in Fig. 1(a), origi- 
nate mainly not from the difference in chemical compo- 
sition, but from that in pnictogen height. By defining the 
contribution t 1 1 nd from the indirect hopping Fe-Pn-Fe as 
the difference between the t\S of the d and dpp mod- 
els (4 nd = tf - t{ pp ), we can estimate f 1 nd to be 0.462 
eV in Ca 4 Al 2 06Fe2P2 and 0.337 eV in Ca4Al 2 6 Fe2 As 2 . 
In our previous work, 24 ) we estimated t 1 1 nd for LaFeAsO 
(ft As = 1-319 A) to be 0.394 eV. These values are reason- 
able because a smaller hp n gives a larger t 1 1 nd for these 
three systems. The direct hopping of LaFeAsO, tf pp = 
—0.243 eV, is much smaller than those of Ca4Al20eFe2P2 
and Ca4Al2 0eFe2As2, which is due to the larger a of 
LaFeAsO. 



Table II. Orbital energies and transfer integrals (in eV) of Fe 
dx 2 — Y 2 an( l P n P Wannier orbitals of dpp models. The descrip- 
tions P and As indicate Ca4Al20eFe2P2 and Ca4Al20eFe2As2, 
respectively, for the chemical composition and crystal structure. 
Transfer integrals obtained from d models are also shown in 
parentheses. The origins of orbital energies £ are set to the re- 
spective Fermi levels. 







Crystal structure 


Composition 




P 


As 




£ X 2 -Y 2 


-0.977(0.023) 


-0.854(-0.204) 




tl 


-0.370(0.092) 


-0.362(-0.044) 




t 2 


-0.080(0.119) 


-0.080(0.065) 




t 3 


-0.014(-0.025) 


-0.011(-0.009) 


P 


£x 


-2.491 


-2.383 




£Z 


-2.537 


-2.846 




tx 


0.153 


0.069 




tz 


-0.646 


-0.537 




£ X 2 -Y 2 


-1.000(0.047) 


-0.882(-0.197) 




tl 


-0.387(0.102) 


-0.374(-0.037) 




t 2 


-0.089(0.121) 


-0.087(0.067) 


As 


t 3 


-0.015(-0.030) 


-0.011(-0.012) 


?X 


-2.496 


-2.391 




£ z 


-2.566 


-2.889 




tx 


0.128 


0.053 




tz 


-0.666 


-0.552 
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3.2 Alloying effect on Ca±Al 2 0§Fe 2 AsP 

To study whether or not even the electronic struc- 
ture of a system in which P and As coexist is de- 
termined mainly by crystal structure, we unfolded the 
band structures and compared them in the three cases 
of Ca 4 Al 2 6 Fe 2 (Asi_ :c P :c )2: Ca 4 Al 2 6 Fe2P2 (x = 1), 
Ca 4 Al 2 06Fe2As2 (x = 0), and Ca 4 Al 2 6 Fe2AsP (x = 
0.5). For the x = 1 and systems, the average crystal 
structure of Ca 4 Al2 0eFe2P2 and Ca 4 Al2 0eFe2As2 was 
used. The fractional height of the pnictogen atoms used 
is thus zp n = 0.0924. For the x = 0.5 system, we used 
a y/2 x y/2 x 1 supercell containing four Fe atoms [see 
Fig. 4(a)]. The lattice parameters and atomic coordinates 
except P and As were set to the same as those for the 
x = 1 and systems. The P and As atoms were set at 
the experimental heights in Ca 4 Al2 0eFe2P2 (^p = 1.306 
A) and Ca 4 Al2 0eFe2As2 (Z^as = 1.500 A), respectively. 
The fractional height of the pnictogen atoms used are 
thus zp = 0.0861 and za s = 0.0989. The dimensions 
y/2 x y/2 x 1 of the supercell are the smallest needed for 
the realization of equal distributions of the two atomic 
species above and below the Fe layer. We unfolded the 
Fe d bands of the x = 0.5 system using eq. (1). Specif- 
ically, we first constructed 20 maximally localized Wan- 
nier functions of the Fe d character and then unfolded the 
interpolated band structure by regarding the supercell as 
consisting of four primitive cells, each of which contains 
one Fe site: a sc = 2a pc ,6 sc = 26 pc , and c sc = c pc . We 
provide in Fig. 4(b) the unfolded band structure for the 
x = 0.5 system as blurred curves and that obtained by 
fixing the weights of the delta functions in eq. (1) con- 
stant as solid curves. It is clearly seen that the unfolded 
band structure is essentially 20-fold and that the weights 
of the delta functions, however, make five of them notice- 
able for almost all the wave vectors. The unfolded band 
dispersion of the x = 0.5 system exhibits different fea- 
tures along the two paths, one from (-zr, 0) to (7r, it) (cor- 
responding to wave vectors along the b direction) and the 
other from (0, tt) to (7r, tt) (corresponding to wave vec- 
tors along the a direction), owing to the inequivalence 
of the a and b directions [see Fig. 4(a)]. Noticeable gap 
openings of three of the unfolded bands are seen on the 
latter path, while the continuity of the unfolded bands 
on the former path looks firmer. They originate from the 
difference between the arrangement of atoms of differ- 
ent species along the a direction, and that of the same 
species along the b direction. 

Figure 4(c) shows the unfolded band structures for 
the x = 0,1, and 0.5 systems. Not only the unfolded 
band structures for the x = and 1 systems, for which 
the exact unfolding is possible, but also that for the 
x = 0.5 system all look quite similar, especially in the 
vicinity of the Fermi levels, as well as in Fig. 3. We also 
confirmed that the band structure of Ca 4 Al20eFe2P2 
(Ca 4 Al20eFe2As2) with an average crystal structure ex- 
cept for the P (As) atom with its experimental height has 
three (two) Fermi surfaces around T. This indicates that 
pnictogen height is a crucial factor for the multiplicity 
of Fermi surfaces around T. It is interesting to see the 
proximity of the Fermi level at (0, 0) for the x = 0.5 sys- 



tem [see Fig. 4(c)]. The top of the second-lowest unfolded 
band sits 0.2 eV below the Fermi level, which is well re- 
produced by the x = and 1 systems with an average 
crystal structure despite the x = 0.5 system containing 
pnictogen height disorder. There is an atomic configura- 
tion for the chemical composition x = 0.5 using the con- 
ventional unit cell (see Fig. 2), in which the atoms of one 
of the two species are all located above the Fe layer and 
those of the other species below the layer. We also cal- 
culated the electronic band structure of this system (not 
shown) with the experimental pnictogen heights of the 
P and As atoms, and found that its unfolded bands near 
the Fermi level are accurately reproduced by the x = 
and 1 systems, as shown in Fig. 4(c). These results cor- 
roborate the independence of the electronic structure of 
a Ca 4 Al206Fe2(Asi_ :c P :c )2 system of the chemical com- 
position. We expect an accurate calculation of a band 
structure near the Fermi level for an arbitrary x using a 
crystal structure modeled simply by linearly connecting 
the Ca 4 Al20eFe2P2 and Ca 4 Al20eFe2As2 crystal struc- 
tures. 

3.3 Generalized susceptibility of 

Ca 4 Al 2 6 Fe 2 (As X - x P x ) 2 
Given the results obtained above, we calculated the 
generalized susceptibility xo of Ca 4 Al2 0eFe2P2 and 
Ca 4 Al20eFe2As2 with gradually varying their crystal 
structures. We expect that the xo calculated in such 
a way accurately reproduces the qualitative behavior 
of Ca 4 Ai206Fe2(Asi_ :c P ;r )2 for arbitrary x. The lat- 
tice parameters and atomic coordinates c used are lin- 
early parametrized by the crystal structure ratio A as 
c(A) = Acp + (1 — A)ca s . A = and 1 correspond 
to the Ca 4 Al 2 06Fe2As 2 and Ca 4 Al 2 06Fe 2 P2 crystal 
structures, respectively. We calculated Xo(q) using eq. 
(2) and plotted it in Fig. 5(a) for Ca 4 Al2 0eFe2As 2 
and Ca 4 Al20eFe2P2 with A = 0,0.5, and 1. Those in 
doped cases with ±0.1 electrons per Fe atom with the 
rigid band approximation are also shown. The plotted 
Xo is seen to be sensitive to the difference in chemi- 
cal composition compared with the band energy. The 
peaks of xo are located at M in the undoped cases of 
both Ca 4 Al20eFe2As2 and Ca 4 Al2 0eFe2P2, while, in the 
doped cases, the peaks are located at incommensurate 
q's close to M. We observe that the heights of the peaks 
are reduced by doping for both Ca 4 Al20eFe2As2 and 
Ca 4 Al20eFe2P2, which suggests that their magnetic in- 
stability is suppressed under doping. Figure 5(b) shows 
the xo(M)'s of Ca 4 Al 2 6 Fe 2 P2 and Ca 4 Al 2 6 Fe 2 As2 as 
functions of A. Although the Ca 4 Al2 0eFe2P2 systems 
with smaller band widths give a larger Xo(M) than the 
Ca 4 Al20eFe2As2 systems, the qualitative behavior of 
monotonic decrease in Xo(M) for both Ca 4 Al20eFe2P2 
and Ca 4 Al2 0eFe2As2 are found to be the same in the 
entire range of A. We found that, with increasing A, the 
number of the a Fermi surfaces around (0, 0) increases 
from one to two near A = 0.75 for both Ca 4 Al20eFe2P2 
and Ca 4 Al2 0eFe2As2 . The larger number of Fermi sur- 
faces for A > 0.75 would be preferable to larger xo be- 
cause of the larger number of scattering channels. The 
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Fig. 4. (Color) (a) Bird's-eye view of an Fe layer of a y/2 x y/2 x 1 
supercell for an x = 0.5 structure is shown on the left. The 
layer viewed along the c-axis is shown on the right, together with 
the primitive lattice vectors of the supercell. Brown, purple, and 
green balls represent Fe, P, and As atoms, respectively, (b) The 
blurred curves are for the unfolded 20 bands of the d model for 
the x = 0.5 system using the supercell. The solid curves rep- 
resent the unfolded bands with weights of delta functions fixed 
constant. The coordinates of high-symmetry points are with re- 
spect to the primitive reciprocal lattice vectors of the primitive 
cell, (c) Unfolded band structures of x = 0, 1, and 0.5 systems. 
The solid and dashed curves are for the 5-band structures of the 
x = 1 and systems, respectively, with the crystal structure 
fixed at the average of Ca4Al20eFe2P2 and Ca4Al20eFe2As2. 
The origins of energy are set to the respective Fermi levels. 



of states. These observations suggest that the effects of a 
small A overcome the difference in the number of a Fermi 
surfaces. 
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calculated xo? however, exhibits the trend of a mono- 
tonic decrease. To examine the factors determining this 
trend, we plot the band energies E xy (7r,7r) in Fig. 5(c) 
for the undoped systems as functions of A. E xy (7r,7r) is 
the height of the d x 2 _y 2 -derived band. As is shown, a 
smaller A gives a higher E xy (ir,ir) and thus a more ef- 
fective 7 Fermi surface, and the difference in chemical 
composition has only a minor effect on the values of 
Exyiji-,^)- Furthermore, a smaller A (a larger a and a 
larger hp n ) leads to narrower bands and a larger density 



Fig. 5. (Color) (a) Generalized susceptibilities Xo(q) of 
Ca4Al2C>6Fe2P2 and Ca4Al20eFe2As2 for crystal structure 
ratios A = 0,0.5 and 1. Those in doped cases are also shown. 
(t>) Xo(M) in doped and undoped cases of Ca4Al20eFe2P2 
and Ca4Al2C>6Fe2As2 as functions of A. (c) Band energies of 
dx 2 —Y 2 -derived bands at (tt,tt) for undoped systems. 



4. Conclusions 

We performed first-principles calculations of 
the electronic structures of Ca4Al20eFe2P2 and 
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Ca 4 Al 2 6 Fe 2 As 2 . Ca 4 Al 2 6 Fe2P2 (Ca 4 Al 2 06Fe 2 As 2 ) 
was found to have three (two) Fermi surfaces around 
T. We found that the systems of the same crystal 
structure but different chemical compositions give rise 
to similar band structures. By constructing the maxi- 
mally localized Wannier functions and analyzing their 
transfer integrals, we demonstrated that the differences 
between the band structures of Ca4Al 2 06Fe 2 P 2 and 
Ca4Al 2 0eFe 2 As 2 originate mainly from the difference in 
pnictogen height. Specifically, the hybridization between 
the pz orbital and its neighboring d x ^-Y 2 orbitals plays 
important roles in determining the electronic structures. 

We analyzed the electronic structure of 
Ca4Al 2 0eFe 2 AsP by unfolding its electronic bands. 
We found that its band structure resembles the band 
structures of Ca4Al 2 06Fe 2 P 2 and Ca4Al 2 0eFe 2 As 2 with 
an average crystal structure, although Ca4Al 2 0eFe 2 AsP 
contains pnictogen height disorder. We therefore ex- 
pected that an accurate calculation of band structures 
for an arbitrary x is possible using a crystal structure 
modeled by linearly connecting the crystal structures of 
Ca 4 Al 2 6 Fe 2 P 2 and Ca 4 Al 2 6 Fe 2 As 2 . 

We then calculated the generalized susceptibility xo 
for Ca4Al 2 06Fe 2 (Asi_ x P a; ) 2 using crystal structures 
constructed from the linear interpolation of the crystal 
structures of Ca4Al 2 0eFe 2 P 2 and Ca4Al 2 0eFe 2 As 2 . It 
was found that a smaller crystal structure ratio A leads 
to a larger xo because of the larger density of states and 
the more effective 7 Fermi surface. 

Further theoretical and experimental investigations of 
the structural and magnetic properties of perovskite-type 
systems are required to clarify and optimize the super- 
conducting properties of such systems. 
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